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Introduction. Loss of stability and buckling of a round plate may be observed if the plate is loaded on the lateral surface. 
The solution to this problem is based on a bifurcation approach. In this case, a plate is considered as a nonlinear elastic 
body. It is important to choose the relation between stresses and deformations in sustainability problems of nonlinear 
elasticity. Simple laws of state (constitutive equations) were considered in early works devoted to this problem, for 
example, material of the “harmonic type” (Sensenig). 

Materials and Methods. Equations of neutral equilibrium for round plates made of Murnaghan and Blatz-Ko materials 
are obtained. Assuming a uniform initial deformation on the plate, the stability problem is considered. Strict three- 
dimensional neutral equilibrium equations provide exploring related forms of equilibrium taking into account physical 
and geometric nonlinearity. Derivation of these equations is based on the application of the theory of superposition of 
small deformation on the final one. 

Results. Progress in solution to the corresponding secular equation (with non-linear parameter entry) for practically 
important laws of elasticity of Murnaghan and Blatz-Ko is possible using the numerical methods only. The developed 
method for calculating bifurcation values of loading parameters makes it possible to analyze the effect of nonlinearity. 
Discussion and Conclusions. The influence of physical and geometric nonlinearity on the upper critical value of the 
initial deformation parameter is explored. The results obtained can be used under the assessment of reliability of elastic 
third-order moduli for various physical materials. Data on these moduli is still scarce. The numerical research has 
shown that the constants given in some sources should be treated with caution. The use of elasticity moduli in the law of 
state of Blatz-Ko is also discussed. 
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Introduction. Currently, the study of new relatively simple and adequate laws of state for various materials 
that require considering nonlinear deformations is an actively developing area of continuum mechanics [1-10]. In the 
framework of the theory of superposition of small deformation on the finite one, three-dimensional equations of neutral 
equilibrium are derived for the materials of Murnaghan and Blatz-Ko. Based on these equations, an example of end 
buckling of a uniformly compressed round plate is considered [11, 12]. 

Materials and Methods 

Equations of neutral equilibrium. Let 7, @, z be the cylindrical coordinates of a point in an undeformed 
state. We assume that the initial deformation of the body is determined by the radius vector R: 

R=R(r)e, + dz1,, (1) 
where R(r)=ar, parameters a, d are constants, e., e,, i, are basis vectors. 

For the coordinates of this point in the initial deformed state, we have: 

R=ar, ®=9, Z=dz. 

Therefore, the strain gradient, the Finger strain measure, and its principal invariants are determined through the 

relations: 
VR =VR’ =a(e,e, + pe) +di,i,, F=VR’, 
Late td, Law sied, Lad ; (2) 

From now on, V, V_ is the nabla operator in the metric of the undeformed and initial-deformed state: 


V=VR"-V. 
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We will use the neutral equilibrium equations proposed by A. I. Lurie [13]: VO=0, where the tensor © is a 
linear differential operator over the vector of additional displacements W. The expressions of the components of this 
tensor concretized with account of the laws of state of Murnaghan and Blatz-Ko were obtained in [14, 15]. In the 
representations of the tensor ©, the components are some functions defined as a result of solving the boundary value 
problem of the initial deformation. 

Research Results. To consider the bending form of the plate equilibrium bifurcation, similarly to [11, 12], we 
accept the following additional displacement vector 

W =u(7,z)e, +W(7, Z)i; - (3) 

Considering (1), (2) and (3), Lurie tensor has the form: 

© = Az,e, + Be,e, + Cit, + Gi,e, + He,i,. 
Here, 
Ou Ow u Ou Ow u Ou Ow 


Aad ee 4 Baa ee ea ee Ba. 
r or Oz r or Oz r or Oz 





G= (eRe g clan H= ee 
Oz or Oz or 


We note that the structure of these operators is typical in the stability problems of cylindrical nonlinear elastic 





bodies. Omitting the expressions of the remaining coefficients, we give, for example, formulas for A,. In the case of 
Murnaghan law of the state [13, 14] A, is expressed through the relation 
4 
AF Oo we | 
d 2 2 
Henceforward, 4 and p are Lame elasticity moduli, v,, v,, v, are the third-order elasticity constants. For 


Blatz-Ko material [15]: 








A= 2u(— B) 20D aye a8 B bam) ; 


aff, 
where f is the refining elasticity modulus. 
The equation of neutral equilibrium is equivalent to the system of differential equations with respect to the 


2: 2 
A,| — wes Pee {4 +— 6, \o "$6, <0, 
or ror or d ordz d 


2 2 
(4410) 2[(Ho4)o1,{ 2 Lae J 12) oY 2G, 
r 


components of the vector W : 





(4) 





dz \ ar or? rr Or) d * 
Assuming that the plate is loaded with uniform pressure along the lateral surface, we supplement the 
differential equations (4) with the equilibrium conditions at the boundary: 
pa, 79 Onl,_,, =O 





At the ends, i.e., at 
h 
z=45 ©,,=0, ©, =0. (5) 


We apply the substitution 


u=X (z)J,(k,r), 6 
w=Z,(z)J,(K,7). 
Here, n=1,2,..., J,, J, are Bessel functions of zero order and first kind, and numbers k,7, are null functions J,(x). 


Let us assume x =k,r and use the equations 


dS (kr) _ d°J,(k,r) 
dr 


; ‘ 7 dJ,(k,r) 
kJ \ (x), ae TO), a 


=k, J(x) = kJ} (x) , 
and well-known identities for Bessel functions 
1 1 1 
fe) =-J,, J," +—J) =-J,- 
x x x 


Having completed the above steps, we obtain that the variables in the differential equations (4) and under the 
boundary conditions (5) are separated. We arrive at the following boundary value problem for ordinary differential 
equations: 
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pa k; 4 Xx, k, © = Le =0, 
H, G, (7) 
zee iz 4p tS yr 9 
3 3 


k,C,X, (£h/2)+C,Z' (+h/2) =0, (8) 

k,G,Z,, (£h/2)-G,X, (+h/2) =0. 
Two homogeneous linear ordinary differential equations (7) and four boundary conditions (8) result in the 
eigenvalue problem with a nonlinear occurrence of the parameter. In this problem, such a parameter is the 


quantity (1 - a) . The system (7) is ended in the standard form: 


Vi =Vo, 
A G,+C 
' 2 2 2 
=k, = tk, OM 
H, G, 
V3 = Vas 
G, G,+C 
y, =k Sy, -k, A 2" 
C, C; 


Te 


The following notations are accepted here: argumentt=z, (3,3.y3¥33)4) =(X,3X75Z,3Z/) . Let the 


fundamental system of solutions be the following four vectors y, = (1,3 ¥3;3 a3 Vai y . Then, the general solution to the 


system is given by (X,,; X/;Z,;Z/ Vode 


4 

=> éy, , where ¢€. are arbitrary constants. For example, we set the initial data 
i=l 

for z=—h/2 by the columns of the fourth-order unity matrix. Then we solve numerically the Cauchy problem with 

these initial conditions. As a result, on the right-hand side, we obtain the values of the basis functions, ie., 

vectors y,(+h/2) . Using the boundary conditions (8), we arrive at a homogeneous system of linear algebraic equations 

with respect to ¢, ¢,, €3, 4: 

4 

La,¢, =9, where i= 1, 2, 3, 4. 

j=l 


Here, the coefficients are the elements of matrix A: 
aj = Kk,Cy,; (—h/2) + C34; (-h/2) » Ag; = KC; (h/2)+ C34; (h/2) 3 
a,, =k,G,y3; (-h/2)- Gyo; (-h/2) » 4,, =k, G,ys; (h/2) -Giyo; (h/2) : 
The homogeneous system has a nontrivial solution if the condition is met 
detA =0. (9) 
The determinant expression includes the load parameters a, d , Bessel null-functions, as well as the elasticity 
moduli 2, uw, v,, v,, v; (for the Murnaghan material) or 2, uw, £ (for the Blatz-Ko material). 


The parameters a, d , which set the initial deformation, are interconnected. The axial force acting on the cross- 
R, 

sectional area is determined by the relation [13]: Q=2z/ o,RdR, where o, is the physical component of the stress 
0 


tensor. Since the initial deformation is assumed in the form (1), the Cauchy stress tensor T and the Finger strain 
measure F are coaxial. 
Moreover, the stress tensor is constant: 
2 fos Bs : ee 
T= ae (a’(e,e, +€,€,) +4 i,i;)-c™ (a*(¢,e, +€,€,) + d*ii,) +c (€e, +e,e, + i,i,)) ; 
The ends of the plate are free of load, thereforeo, =0. This implies the condition relating the coefficients of 
the Finger law 
co 
c —c%d? + =0. (10) 
d 
Considering the law of state, we can write a specific expression of the condition (10), which establishes a 
connection between a and d. So, for the Murnaghan material, we get the condition: 


[Sis Bay, Ja! +[a42u-Shsy, —4y, +(y, anya \a - 
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3 15 
a 4 AOS ia eS. 
4 2 4 8 
In particular, if vy, =v, =v, =0, then we get 
Pie (1-7). 
+2 
When considering the Blatz-Ko law of state, d is expressed by the formula 


dea. 
Thus, through setting the elasticity moduli and Bessel null-functions, we find the bifurcation values of the 





initial strain parameter a, from (9). 
Note that the value A =1-a, is a small parameter for relatively thin disks. So, in the classical theory of plate 
buckling, the critical value of a, is determined by the formula 
(3,8317)7 (4) 
a, =1-— ; 
° 12(1+v) t ) 


where v is the Poisson's ratio, and 3.8317 is the first root of the Bessel function J,(x). 





When solving the initial boundary value problem, it is assumed that the Signorini’s perturbation method can be 
applied. At this, we assume that the coefficients of the operator © depend on a small parameter A in a power-law 
manner. This means that the boundary conditions for the incremental displacement components A are specified by the 
matrix A = A(A). In this problem, a partial eigenvalue problem — the determination of the lowest eigenvalues — has 
physical meaning. The higher degrees A , which contain the matrix elements, slightly affect the values of the smallest 
roots of the secular equation (9). If we restrict ourselves to the linear theory under solving the initial problem, then 
A(A) is a regular binomial [16]. 

Numerical experiments under studying the stability of nonlinear elastic bodies of not too large relative 
thickness confirm this conclusion. Therefore, a characteristic equation with non-linear occurrence of parameter (9) can 
be replaced by the characteristic equation of the linear operator. Iterative processes that converge to one eigenvalue, 
where a number close to the value in the theory of plates is chosen as the null approximation, can be applied. 

Discussion and Conclusions. Using the equations obtained above, a numerical analysis of the influence of 
physical and geometric nonlinearity on the value of the upper critical parameter a, is performed. The calculations are 
implemented in the Matlab environment for various options of specifying elasticity moduli, relative plate thickness, and 
waveformation number 7 . It is found that in all cases of loss of disk stability, the first-mode buckling that corresponds 
to the minimum critical value of the parameter, which corresponds to the root of the Bessel function 3.8317 [17], 


answers the smallest critical parameter value a, . 
Table 1 shows the critical parameter values a, for plates with a relative thickness of 0.05; 0.1; 0.15; 0.2; 0.25; 
and 0.3. 























Table 1 
Critical parameter values a, for plates of various relative thicknesses 
ne Relative plate thickness 
" 0.05 0.1 0.15 0.2 0.25 0.3 
I 0.9922 0.9869 0.9800 0.9660 0.9530 0.9330 
0.9972 0.9867 0.9735 0.9549 0.9300 0.8914 
> 0.9977 0.9908 0.9801 0.9664 0.9504 0.9332 
0.9967 0.9869 0.9690 0.9410 — — 
3 0.9976 0.9907 0.9800 0.9662 0.9502 0.9330 
0.9965 0.9856 0.9665 0.9372 — — 
x 0.9977 0.9907 0.9798 0.9659 0.9497 0.9324 
0.9881 — — 0.9965 0.9613 0.9861 
5 0.9914 0.9865 0.9783 0.9668 0.9520 0.9350 
0.9929 0.9729 0.9341 0.9114 — — 
6 0.9978 0.9912 0.9800 0.9675 0.9520 0.9340 
0.9985 0.9941 0.9871 0.9775 0.9651 0.9539 
7 0.9985 0.9941 0.9859 0.9735 0.9560 0.9320 
0.9985 0.9941 0.9863 0.9749 0.9594 0.9396 
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The following numbers indicate the materials listed below. 
Steel Rex 535 (A =1.09, w=0.818, v, =-1.75, v, =-2.40, v, =—1.69). 


Steel SOHGSM2F (24=1.129, w=0.803, v, =—-2.8, v, =—-2.7, v, =—1.87). 
Steel Hecla 37 (A=1.11, w=0.821, v,=—3.58, v, =—2.82, v,=~-1.77). 
Steel Hecla ATV (A=0.87, w=0.716, v,=0.34, v, =—-5.52, v,=—-1.0). 
Beryllium bronze (A=1.042, u=0.49, v,=—-4.0, v,=-1.7, v, = -0.6). 
Organic glass (A=0.39, uw=0.186, v,=—0.078, v, =-0.07, v, =0.047 ). 


In the first six cases, the Murnaghan law of state is considered. The last option presents the results for the 
material of Blatz-Ko. The top number in the table cell refers to the case in which physical non-linearity is not taken into 


ON Hs Be ke 


account, Le., v, =v, =v, =0 in the Murnaghan law; £ =0 in the Blatz-Ko law. The second (lower) number takes into 
account physical nonlinearity. In the Murnaghan law, data were used for the third-order elasticity moduli from [14] in 


units 10’ ale 
cm 
In the last row of Table 1, there are the results for the Blatz-Ko law, which was chosen in a simplified version 


(the Knowles-Sternberg equation): the Poisson's ratio was taken equal to 0.25, and the refining module f =0.45. A 
dash means no critical value a, was found. 


The analysis of the results provides drawing some conclusions. At small relative thicknesses of the disk, the 
exact theory and the linear theory of plates give the same critical parameter values a,. Geometric nonlinearity has 


marked impact at relative thicknesses greater than 0.1. Physical nonlinearity is even more pronounced. However, care 
should be exercised in choosing third-order elasticity moduli. For example, in the fourth version (Hecla ATV steel) and 
for thin plates, no reliable critical values of the parameter a, were found; although acceptable values indicating a loss of 


stability are observed for the same Lame elasticity moduli. As in the stability problem for a nonlinear elastic sphere 
made of Blatz-Ko material [15], the constant B weakly affects the critical value a, . 
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